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functions, but carried out on a simple grid at a low computational cost, comparable 
to the inexpensive predictive techniques. 
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1. Introduction 

Advanced LIGO gravitational-wave (GW) detectors are expected to start science 
observations in 2015, with Advanced Virgo following shortly after [1, 2, 3, 4, 5]. A prime 
source predicted to be observable by these detectors are compact binary coalescences 
(GBGs) during which neutron star (NS) or black hole (BH) components of the binary 
are driven by GW emission through an inspiral into a merger. Binary neutron stars 
(BNS) have already been observed as pulsars in Galactic binary systems evolving 
due to emission of GWs [6, 7]; extrapolations from these observations yield estimated 
coalescence rates of BNS systems as detectable by the Advanced LIGO-Virgo network 
([8] and references therein). 

1.1. Binary coalescence model 

The information available in an observed circular GBG event is encoded through the 
waveform h, fully described by a 15-dimensional parameter vector 9 which in turn can 
further be divided into two parameter sets. The amplitude and phase evolution of 
the waveform are governed by the intrinsic parameters which include the masses and 
spins of the binary components. The extrinsic parameters describe the projection of the 
gravitational wave emitted by the binary onto the detectors and include the object’s 
position and orientation on the sky, its luminosity distance di as well as the time and 
phase of the waveform at coalescence {tc and ff). The level of correlation between the 
two parameter sets is only marginal, apart for tc and fc which are strongly correlated 
with intrinsic parameters. Therefore we will simplify the analysis by hxing the remaining 
extrinsic parameters at their true values without signihcant impact on the recovery of 
the intrinsic parameters; we return to this point in section 4. 

In this study we consider as an example a BNS with component masses mi = 
1.45Mq and m 2 = 1.35Mq. While there is a strong degeneracy between the mass and 
spin parameters, which affects their measurement uncertainty [9, 10, 11, 12, 13], the 
BNS components are in this study taken as non-spinning for demonstration purposes, 
as our focus here is on parameter estimation techniques rather than astrophysical 
predictions. Therefore, the original 15-dimensional parameter space is reduced to only 
four parameters, the two masses as well as tc and 0c; without loss of generality in the 
method applied. The waveform approximant used is TaylorF2 [14] which for a face-on 
binary directly overhead the detector describes the frequency-domain waveform as 

assuming no higher-order corrections to the waveform amplitude. Here, Aic = 
{mim 2 )^^^{mi m 2 )~^^^ is the chirp mass and the phase 4/(/, 9) is given as 

'!'(/, ^) = 27r/G - 0c - ^ + + Pk log(z/))z/'' 

' k=0 


^7r-2/3M5/6^-7/6g*vI.(//)^ 


( 1 ) 


( 2 ) 
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where u = (7r(mi + m 2 )f Y^^ and the coefficients ak, Pk are functions of Aic, the spins 
(here assumed to be zero) and the mass ratio p = {mim 2 ){rni + 7712 )“^ for each post- 
Newtonian order /c/2. 

The likelihood of observing data d = hr + n, containing the true waveform hr and 
stationary, Gaussian instrumental noise n, given a waveform model h{6) described by 
parameters 6 is 

L{9) oc exp ^“2 + hr — h{9)\n -\- hr — h{9))^ . (3) 

The term in the angle brackets is the noise-weighted power in the residuals after 
subtracting the assumed model from the data. The angle brackets denote an inner 
product between the Fourier transformed time series d{f) and the complex conjugate 
of b{f) weighted by the noise power spectral density 5'„(/): 


{a\b) = 4 Re 


Snif) 


df. 


( 4 ) 


The likelihood (3), together with the specification for the waveform model and the 
noise power spectral density 5'„(/), which is chosen to match the Advanced LIGO zero¬ 
detuning high-power design sensitivity [15], constitutes the model H. 

For this study we will consider a BNS observed directly overhead in one detector 
at a signal to noise ratio p = 12, defined as in [16]: 

{d\h{9)) 


p = max 
e 


{h{9)\h{9)) 


( 5 ) 


The injected noise realisation in our example is chosen to be n = 0, which hxes the 
maximum likelihood waveform to be Ht- 


1.2. Bayesian inference 


Given the model H as described in subsection 1.1, (3) can be viewed as a probability 
density function (PDF) p{d\9,H) for observing data d = hr + n given parameters 9. 
Together with a PDF on the parameters detailing the prior knowledge and expectation 
about 9 under FT, p{9\H), a posterior PDF can be obtained through Bayes’ theorem as 


p{9\H)p{d\9,H) 
p{d\H) 

The term in the denominator is called the evidence and can be used for comparing 
the effectiveness of different models’ ability to describe the data; for this study, it is 
a normalizing constant for the posterior PDF. The task of parameter estimation can 
therefore be dehned as hnding accurate and robust methods for evaluating the posterior 
probability given by (6) in order to obtain the accuracy with which different parameters 
can be recovered. This accuracy is usually quantihed in terms of credible regions GRp 


p{9\d, 
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containing a fraction p of the total posterior probability. These credible regions are 
dehned as 

CRp = min A (7a) 

such that the parameter-space area A satishes 

p = [ p{9\d,H)d9. {7b) 

J A 

The posterior PDF will have the same dimensionality as 9 and a shape strongly 
dependent on the complexity of H. For the PDFs encountered in CBC data analysis, 
very intricate structures with high levels of multimodality and parameter degeneracies 
are common, and their prevalence has spurred the development and use of advanced 
analysis tools for sampling the posterior PDF. 

1.3. Stochastic sampling 

As part of the publicly available LSC Algorithm Library (LAL) [17], the Bayesian 
framework in LALInference implements several methods to stochastically traverse the 
parameter space [18, 19], creating a set of individual samples 6*j distributed according 
to the posterior PDF in (6). For this study, we used a parallel-tempered Markov 
chain Monte Carlo (MCMC) algorithm [20, 21], implemented in LALInference as 
LALInf erenceMCMC. The number of MCMC samples required to fully describe the 
posterior PDF scales only weakly with the number of dimensions, allowing for efficient 
exploration of high-dimensional parameter spaces. However, only a small fraction of the 
collected samples will be statistically independent, due to imperfect jump proposals in 
the complex multimodal parameter spaces. In order to accurately represent the posterior 
PDF for the 4-dimensional parameter vector 9 from subsection 1.1, (P(IO^) independent 
samples are necessary [22]. This requires a total of (P(10®) LALInf erenceMCMC posterior 
samples, and therefore 0{19S) individual likelihood calculations spread across the 9 
differently tempered parallel chains used in this example. 

l.f. Paper organisation 

In section 2 we present the “cumulative marginalized posterior” method for accurately 
and efficiently evaluating credible region contours on a grid. We show that credible 
regions can be computed with a suitable choice of a low-density grid of samples, 
described in subsection 2.2. This makes it possible to produce credible regions with 
the same accuracy as given by the stochastic samplers, but at a greatly reduced 
computational cost, as shown in subsection 2.3. In section 3 we discuss alternative 
grid-based methods for estimating credible regions, including the iso-match contour 
technique advocated by Baird et al. [23], and show their relative performance against 
the cumulative marginalized posterior method (subsection 3.3). We conclude in section 4 
with a discussion of these results and a proposal for implementing the cumulative 
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marginalized posterior method in a low-latency parameter estimation framework for 
GW signals from CBCs. 

2. Discretizing the credible regions 

As an alternative to the stochastic sampling methods we propose a method to evaluate 
the integral in (76) using samples at pre-determined coordinates in parameter space. The 
simplest implementation places samples into centres of pixels distributed in a uniform 
rectangular grid and approximates the posterior as constant within a pixel. The choice 
of the sampling grid is critical to the efficiency of this method, and will be discussed in 
subsection 2.2. 

The required number of likelihood calculations, and therefore the required number 
of pixels, scales exponentially with the number of dimensions in the grid, so in order to 
minimize the computational requirements a modihed version of (3) was implemented. 
We use a likelihood function marginalized over tc and fc, thus removing the need for 
sampling those dimensions without affecting the PDF in the mass parameters of interest. 
This marginalized likelihood function is presented in [24], which also provides constraints 
on the resolution at which the waveforms must be sampled in both the frequency and 
time domains, in order to guarantee no loss of information when computing the inner 
products. 

2.1. Cumulative posterior on a grid 

Taking advantage of the natural [Adc, p\ parameterization of the waveforms, we 
constructed a uniformly spaced rectangular grid across these parameters. For each 
pixel in the grid, the likelihood corresponding to the value of 9 at its midpoint was 
calculated. We assumed a prior PDF on the mass distribution of BNS systems to be 
uniform in [mi, m 2 ] rather than [Aic^rj] (transformed appropriately when working in 
the [Adc,"/]] space). The product of the likelihood and prior yield the numerator of 
(6). The posterior PDF is obtained by normalizing this quantity by the evidence, the 
denominator in (6), which is approximated as the sum of the likelihood-prior products 
over all pixels in the grid. The posterior PDF is assumed to be valid not just locally at 
9 but instead across a whole pixel. Credible regions can then be dehned by the set of 
pixels containing a fraction p of the posterior PDF, accumulated when traversing the 
pixels in order of decreasing posterior values. 

2.2. Grid placement 

To minimize computational cost and to ensure an accurate representation of the credible 
regions dehned by the integral in (76), the grid samples must be placed as sparsely as 
possible, subject to two constraints: (i) a sufficient fraction of the parameter space 
with signihcant posterior support is covered to enable accurate normalization of the 
posterior; and (ii) the error introduced by approximating the prior-likelihood integral 
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over any pixel as the product of the prior and likelihood at the centre and the pixel area 
is within desired bounds. 

The required extent of the grid can be quantihed in terms of the Mahalanobis 
distance r{9) dehned as 

r = f (g- (8) 

for a set of pixel coordinates 6 spanning a multivariate iV-dimensional PDF f{6) = 
p{6\d,H) with mean fl and covariance S [25]. When f{9) is a bivariate Gaussian, the 
associated cumulative density function is given by <h(r) = 1 — Hence, for a 

maximum error e = 1 — <h(r) in the evidence contained within the grid, the grid must 
minimally contain the pixels bounded by a distance 

Tfe = V—2 Ine (9) 


away from the maximum (mean) of the PDF. The main purpose of this analysis is to 
construct credible regions whose p-value is known to an accuracy no worse than that of 
a stochastic sampler, which is of order 1% for the (P(IO^) samples we typically have (see 
subsection 2.3). We therefore set e = 0.005, and will correspondingly cover the region 
Tfe < 3.25 with a grid. 

The minimum density of pixels within this bound is set by requiring that the 
approximate PDF, computed by discretely evaluating the posterior on a grid, is a 
sufficiently good approximation to the integral {7b). Here we use a very simple 
approximation, namely, we evaluate the integral as a Riemann sum over equal-sized 
pixels, setting the contribution of each pixel to the integral equal to the product of the 
pixel area and the value of the PDF at the centre of the pixel. In this case, a sufficient 
- but not necessary - condition for the total error on the PDF integral to be bounded 
by e is for the fractional error in each pixel to be smaller than e. For an iV-dimensional 
PDF f{9), this fractional difference across a pixel centred at 9q is 




( 10 ) 


where A is the pixel width. For a one-dimensional Gaussian distribution f{9) with zero 
mean and variance the integral in (10) can be represented as a Taylor series: 

/‘0O+A/2 

/ f{0)d9 

Jeo-A/2 

r‘eo+A/2 / n 1 n2 

- ^o) - gp-AO - 9o)^ + ^A9 - 9,)^ + 0{{9 


f{0o) 1 


I do-A/2 




2ct2 


2a4 


9^)^' 


d0(ll) 


/(^o 


A- 


A3 

24^2 


24a^ 


where the hrst non-zero correction term enters at the second order of the Taylor series 
since the approximated PDF is evaluated at the centre of the pixel. The integral will 
be dominated by the last term in (11) for increasing I^qI; hence, the most stringent 
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requirement on the pixel size will come from 6o at the bounds of the integration domain. 
As discussed above, our analysis is restricted to 0 < \0Q\/a < r?,, so (10) becomes 


1 ) 


24 




< e, 


( 12 ) 


Hence, for e = 0.005, the uniform grid size is A ^ O.la, i.e., a total of ~ 60 pixels are 
required per parameter-space dimension. Therefore, assuming no correlations between 
parameters, a grid of ~ 3500 pixels is required to achieve 99.5% coverage of the posterior 
region and sub-percent net credible region identihcation. 

In practice, the grid size can be signihcantly reduced in a number of ways. The grid 
size need not be regular; rather than requiring a hxed fractional error per pixel, we could 
require a hxed contribution to the absolute error, which would allow us to increase the 
size of pixels near the bounds of the integration volume that contain a very small fraction 
of the PDF but set the most stringent requirements if the fractional error criterion 
is used. More accurate integration can be obtained by higher-order schemes, such as 
Simpson’s rule, lowering the minimum number of grid points necessary to achieve a given 
accuracy. The grid need not be rectangular, but could be an iV-dimensional sphere of 
dimensionless radius rh, achieving a signihcant volume reduction in a high-dimensional 
space over a cube enclosing such a sphere, as assumed above. 

If the assumption of uncorrelated parameters is relaxed, (9) and (12) will still 
give the number of pixels required, but their relative placement needs to be altered. 
A misalignment between the grid and the PDF caused by correlated parameters or a 
non-ellipsoidal posterior PDF will reduce the validity of the previously given scaling 
relations. This can be overcome by either oversampling the grid; a coordinate rotation 
to align the grid and the PDF; or a dynamical placement of the pixels, adapting the local 
pixel density to a preliminary PDF estimated from a coarse grid across the parameter 
space. These extensions to the analysis will be investigated further in future work. 


2.3. Key results 

For the example BNS system described in subsection 1.1, the grid-based cumulative 
marginalized posterior calculation can determine the credible regions of the posterior 
PDF to the same accuracy as LALInf erenceMCMC using only a small fraction (~ 0.1%) of 
the computational cost of the stochastic sampler. This is illustrated in Figure 1, which 
shows 2945 independent samples of the posterior PDF produced by LALInf erenceMCMC, 
colour coded by the credible region they fall into as given by the grid-based cumulative 
marginalized posterior. The credible regions are compared to the fraction of MCMC 
samples falling into the pixels contained within them, i.e., to the credible regions as 
estimated by the stochastic method. The number of MCMC samples falling within a 
claimed credible region is governed by a binomial distribution such that the uncertainty 
in the fraction of samples in CRp is a/p(1 — p)/N ; e.g., for N = 2945 and p = 0.3, the 
uncertainty is 0.8% - consistent with the observed fluctuations in Figure 1. The credible 
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Figure 1. MCMC samples (dots) in [McV] space are colour-coded by the cumulative 
marginalized posterior credible region they fall into. The legend compares these 
credible regions against the fraction of MCMC samples falling into them, which 
corresponds to the stochastic estimate of the fraction of posterior contained within. 
The true parameters of the evaluated BNS system are shown at the turquoise cross 
(corresponding to mi = 1.45Mq and m 2 = I.SSMq). 


regions, and associated uncertainties, estimated for all p G [0,1] are shown in Figure 4 
as a complement to the discrete set of credible regions displayed here. 

3. Comparison with alternative methods: which approximations are 
warranted? 

We have demonstrated that the cumulative marginalized posterior method is both 
accurate and computationally efficient with respect to stochastic samplers. We now 
explore which additional approximations can be made to simplify the analysis further; 
in the process, we investigate the validity of approximate techniques proposed by Baird 
et al. and Hannam et al. [23, 13]. 
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Figure 2. MCMC samples (dots) in [McV] space are colour-coded by the cumulative 
maximized likelihood credible region they fall into. The legend compares these 
credible regions against the fraction of MCMC samples falling into them, which 
corresponds to the stochastic estimate of the fraction of posterior contained within. 
The true parameters of the evaluated BNS system are shown at the turquoise cross 
(corresponding to mi = 1.45Mq and m 2 = I.SSMq). 


3.1. Cumulative likelihood 

In the Bayesian formalism used here, it is often assumed that the majority of the 
information about the posterior PDF originates from the likelihood function alone, 
with only a weak dependance on the prior PDF. To test this assumption, we repeated 
the analysis performed in section 2, but without the inclusion of the prior detailed 
in subsection 2.1. Evaluating the cumulative marginalized likelihood across the same 
[A^c, h] provided only negligible computational savings compared to the cumulative 
posterior as the overwhelming fraction of the computational cost is due to the likelihood 
calculations. The removal of the prior radically changed the shape of the credible regions 
from what was observed for the cumulative marginalized posterior in Figure 1 to the 
more ellipsoidal features shown in Figure 2; moreover, the credible regions computed 
via the cumulative marginalized likelihood method no longer match the posterior PDF 
as evaluated with the MCMC sampler. 
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We next evaluated the effect of marginalization over tc and fc by replacing the 
previously used likelihood function with one which instead maximizes over tc and 
^c, for both the cumulative posterior and likelihood methods. This carries greater 
computational savings compared to ignoring the prior, as maximizing the likelihood 
function reduces the computational time by ~ 1/3 compared to marginalizing over 
tc and fc- However, replacing the correct marginalization with the computationally 
cheaper maximization produces credible regions which are no longer consistent with the 
posterior PDF estimated from MCMC calculations within the uncertainty discussed in 
subsection 2.3 (cf. Figure 4). 

While these two simplifications, particularly the use of maximization in lieu of 
marginalization, can yield reductions in computational complexity, the discrepancies 
introduced with respect to the credible regions produced by either LALInf erenceMCMC 
or the cumulative marginalized posterior, are found to be outside the required tolerance 

level. 


3.2. Iso-Match contours and the Linear Signal Approximation 


As an alternative approach for estimating credible regions and predicting parameter 
accuracy, Baird et al. [23] introduced a method, later implemented by Hannam et al. 
[13], based on the iso-match contours. This method relies on the validity of the Linear 
Signal Approximation (LSA) [26, 27]. Under the LSA, waveforms are assumed to vary 
linearly with parameters, allowing a first-order expansion 

h{e) = hT + hiA9 ^, (13) 


where hi is the partial derivative of the waveform with respect to the parameter and 
A0* = — 9f. When combined with (3), this yields the likelihood function 


L{9) (X exp ( —-{hi\hj)A9^A9^ 


(14) 


assuming n = 0, expressed as a multivariate Gaussian centred at the true parameters 
9t with covariance matrix {hi\hj)~^. 

The method of Baird et al. uses the waveform match M between waveforms Lt 
and h{9) defined as 


M 


max 


{hTvm 

f{hT\hT){hmh(g)) ’ 


(15) 


again maximizing over tc and 0c, as a proxy for the likelihood function. By assuming 
that the LSA is valid, Baird et al. approximated credible region boundaries as contours 
of constant match. 


Mp 


X^(l -p) 

2p2 


(16) 


via an iV-dimensional distribution where N is again the number of dimensions of the 
parameter space remaining after maximization. In addition to using maximization in 
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Figure 3. MCMC samples (dots) in [Mc,g] space are colour-coded by the credible 
region they fall into as determined by the iso-match contour bounding them. The 
legend compares these credible regions against the fraction of MCMC samples falling 
into them, which corresponds to the stochastic estimate of the fraction of posterior 
contained within. The true parameters of the evaluated BNS system are shown at the 
turquoise cross (corresponding to mi = 1.45M0 and = l.SSM©). 


lieu of marginalization, the validity of this approximation relies on the Gaussianity of 
the posterior, and does not include a priori information. 

Calculating the matches given for each pixel in the same [Aic, v] grid as used in 
Figure 1 and 2, we dehned credible regions as the pixels bounded by the iso-match 
contour in (16). Figure 3 compares credible regions given by iso-match contours 
against estimates from the fraction of MCMC samples falling within those contours; 
the differences between the two are statistically signihcant. 

While the matches used by this method are calculated exactly, the restrictions 
implied by the LSA will lead the method to fail if the posterior PDF under investigation 
exhibits even a moderate level of non-Gaussianity. This can originate in the likelihood 
itself, or from the neglected contribution of the prior. This becomes clear for the BNS 
system evaluated here from the high degree of similarity between Figures 2 and 3, 
indicating the validity of the LSA for this system. In this particular case, the Gaussianity 
of the likelihood in [Aic, p] space means that the posterior would have been Gaussian 
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if the priors were flat in [Aic, v] space, so the method could have performed relatively 
well; it does not perform well for flat priors in [mi, m 2 ] space, as indicated by Figure 4 
(see below), because the posterior in this case is no longer Gaussian. 

3.3. Comparasion 

We compare all of the grid-based methods described above in Figure 4. We show the 
differences between the fraction of MCMC samples contained within the various methods 
for estimating credible regions corresponding to credible level p, and the value of p, for 
p G [0,1]. Perfect agreement would correspond to a horizontal line at a deviation of 
zero. However, the hnite number of stochastic MCMC samples from LALInf erenceMCMC 
leads to statistical fluctuations in the deviation; their expected magnitude is indicated 
by an error ellipse (see subsection 2.3) corresponding to one-a fluctuations. 

The cumulative posterior method, using a likelihood function marginalized over tc 
and (l)c, successfully estimates credible regions (apparent deviations at p < 0.3 could 
be statistical, or may be due to the need for sub-pixel resolution to resolve such small 
credible regions). The comparison also further solidihes both the validity of the LSA 
for this system and the effect of the prior on the ability to recover consistent credible 
regions with respect to LALInf erenceMCMC. 

4. Conclusions and future directions 

We have evaluated several grid-based methods for approximating the parameter credible 
regions for a CBC event, and compared these to regions estimated by the stochastic 
sampler LALInf erenceMCMC. We found that evaluating the cumulative posterior on a 
relatively low-density grid allowed us to estimate credible regions to within the statistical 
uncertainty of the stochastic sampler at a small fraction of the computational cost 
(~ 0.1%), while marginalizing over the time and phase parameters and incorporating 
an arbitrary prior. 

On the other hand, ignoring the prior or maximizing over tc and fc instead of 
marginalizing over them introduces a discrepancy in the recovered credible regions with 
respect to LALInf erenceMCMC. In addition, we have demonstrated that the iso-match 
method proposed by Baird et al. is overestimating the credible regions in [A4c, v] space 
compared to a full Bayesian analysis. 

The analysis has been performed on a binary observed in one detector at a 
hxed overhead and optimally oriented position and at a hxed distance giving p = 
12. Compared to an analysis of the same binary using a three-detector observation 
comprising data from the two LIGO observatories and the Virgo observatory (all 
operating at an the same sensitivity as assumed in subsection 1.1), where extrinsic 
parameters describing sky location, inclination, orientation and distance are allowed to 
vary, the recovered two-dimensional credible region in [A4c, rj] space, is not signihcantly 
altered (see Figure 5). The analysis is implemented with the same [AiciV] grid 
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Figure 4. Difference between the fraction of MCMC samples from the posterior PDF 
falling into credible regions predicted by the various methods described above and the 
expected credible level p, as a function of p. Continuous relations for p G [0,1] of the 
data in in Figure 1, 2 and 3 are shown by five curves, corresponding, from the top down, 
to the cumulative marginalized likelihood, cumulative maximized likelihood, iso-match 
contours, cumulative marginalized posterior, and cumulative maximized posterior. The 
filled ellipse indicates the expected 68% level of fluctuation in the fraction of MCMC 
samples due to the finite number of samples. 


as in previous figures, a grid designed for a two-dimensional analysis as described 
in subsection 2.2. The discrepancy, if any, introduced by allowing for eventual 
correlation caused by the inclusion of extrinsic parameters is found to be within the 
statistical uncertainty from the limited number of samples from LALInf erenceMCMC 
(c.f. subsection 2.3). 

The cumulative marginalized posterior method can easily be extended to include 
other parameters, especially spin [13], but as the computational cost scales exponentially 
with the number of parameters, the cost quickly approaches and exceeds the 
computational requirements of the stochastic sampler implemented in LALInf erence. 
Our simple grid-based sampling implementation would be computationally competitive 
with the stochastic sampling for parameter spaces with up to four non-marginalized 
dimensions. However, the computational cost of the grid-based sampler could be further 
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Figure 5. MCMC samples (dots), now drawn from a 9-dimensional analysis, in 
[A4c, g] space are colour-coded by the cumulative marginalized posterior credible region 
they fall into, using the same two-dimensional grid as in previous figures. The 
legend compares these credible regions against the fraction of MCMC samples falling 
into them, which corresponds to the stochastic estimate of the fraction of posterior 
contained within, now allowing for effects of correlation against extrinsic parameters. 
The true parameters of the evaluated BNS system are shown at the turquoise cross 
(corresponding to mi = 1.45M0 and m 2 = 1.35Mq). 


reduced through more efficient grid placement and more accurate integration algorithms. 
Moreover, while standard stochastic samplers such as LALInferenceMCMC are serial 
processes, all pixels in the sampling grid are completely independent, therefore trivially 
allowing for massive parallelization of the cumulative marginalized posterior analysis. 

Through these properties we envision the grid-based sampling method, using a 
cumulative marginalized posterior, to be implemented as a low-latency parameter 
estimation tool for the intrinsic parameters of a CBC candidate event, similar to the 
implementation of bayestar for the extrinsic parameters [10, 28]. In practice, we 
won’t know the true signal parameters which are needed for efficient and accurate grid 
placement. However, we can take the parameters of the highest-match template from 
the search pipelines used for the detection of CBC events [29, 30], which by design of the 
template banks will generally have M > 0.97 [31], as the central point of the grid. The 
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size of the grid can be initially estimated by comparing the SNRs reported by adjacent 
templates in the template bank. Subsequently, the grid can be adaptively rehned: while 
we used a grid with uniformly spaced pixels for this study, the computational cost could 
be reduced further by implementing a non-uniform grid. The reduction in the absolute 
number of pixels required within the grid can enable the inclusion of additional non- 
marginalized dimensions while retaining the computational competitiveness against the 
stochastic sampling methods. 

Moreover, predictive methods such as the iso-match method [23, 13] or the effective 
Fisher matrix approach [32, 33] can be used in combination with the grid-based 
cumulative marginalized posterior technique to provide PDF estimates for dynamically 
laying out the grid. The cumulative marginalized posterior could also be implemented 
as a jump proposal for the intrinsic parameters as part of the stochastic samplers in 
LALInf erence. These possibilities will be explored further in future work. 

In addition, it is important to note that even though the method of parameter 
estimation with the grid-based sampling using a cumulative marginalized posterior has 
been presented here in the context of gravitational-wave astrophysics, the method itself 
is completely general and can be effective whenever the dimensionality of the parameter 
space is sufficiently small to make it competitive with stochastic samplers. 
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